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Abstract 
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1 Introduction 



Understanding fermion masses and mixing remains one of the most important outstanding 
problems in particle physics. In particular, the issue of possible neutrino masses and associ- 
ated lepton mixing is of fundamental interest. Although there is no definite direct evidence 
for nonzero neutrino masses^ they are expected on general grounds: given only the known 
left-handed neutrino fields and the usual Higgs field(s) in the standard model (and super- 
symmetric extensions thereof which stabilize the hierarchy), nonzero neutrino masses result 
generically from higher-dimension operators which are expected to OCClir clt cl SCcllc near to 
that of quantum gravity|, suppressed by associated inverse powers of the (reduced) Planck 



mass, Mp = y he/ (8tiGn) = 2.4 x 10 18 GeV. For example, gauge-invariant dimension-5 oper- 
ators of this type could produce neutrino masses of order m u ~ av 2 /Mp, where v = 250 GeV 
is the scale of electroweak symmetry breaking, and a is a dimensionless constant. (Here Mp 
is an approximate upper bound on the mass which suppresses such operators; it is possible 
that new physics occurs at some intermediate mass scale v << Mj < Mp in such a way that 
dimension 5 operators of this type would give rise to neutrino masses of order a'v 2 /Mj.) Thus 
one can understand on general grounds why neutrino masses are so small. It is not known 
whether there exist any electroweak-singlet neutrino fields. If they do exist, then they could 
lead, via renormalizable, dimension-4 operators, to neutrino masses m u ~ v 2 /Mp, where the 
scale Mp of the electroweak-singlet neutrino mass is naturally » v, again yielding, albeit 
for a different reason, very small m v ||. 

In turn, a natural concomitant of (nondegenerate) neutrino masses is lepton mixing, 
which is thus also a generic expectation.^ It is thus of interest to study models for the 
leptonic mass matrices in the charge Q = and Q = — 1 sectors. The diagonalization of 

1 For reviews and current limits, see Refs. |jj and ||. The apparent solar neutrino deficit is the most 
suggestive indirect evidence at present. The situation with atmospheric neutrinos is unclear. 

2 This expectation is based on the general consensus that pointlike theories of quantum gravity (in par- 
ticular, supergravity) are non-renormalizable and is borne out by the E << (a') -1 / 2 limit of string theories. 

3 The lepton mixing angles are functions of ratios of elements of neutrino matrix elements and of charged 
lepton mass matrix elements, and even though left-handed neutrino masses are small, some of these ratios 
could, in principle, be 0(1). However, a set of conditions for natural suppression of observable lepton flavor 
violation were formulated, and it was shown that the standard model (generalized to include nonzero m„) 
satisfies these j|. 
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these matrices determines both the masses and the observable lepton mixing. In addition to 
the search for neutrino masses and mixing in solar and atmospheric neutrino experiments, 
searches continue at accelerators. In particular, the CHORUS and NOMAD experiments at 
CERN are currently looking for — > v T oscillations[]|| 

In analyzing models of fermion masses and mixing, an important step is to determine 
the number of real amplitudes and unremovable, and hence physically meaningful, phases 
in the mass matrices. Recently, we reported a general solution to this problem J7|. Here 
we extend our explicit construction of rephasing invariants and give further applications to 
various models. We also compare the situation concerning unremovable phases and associ- 
ated invariants in the leptonic sector with that in the quark sector, for which we have also 
recently presented a general analysis ||. 

The organization of this paper is as follows. In section 2 we briefly review our theorem on 
the number of unremovable phases. In section 3 we deal with the question of which elements 
of the mass matrices can be made real by rephasings. We discuss the connection between 
this and the complex rephasing invariants. This section includes explicit expressions for the 
sixth order complex invariants. In section 4 we apply our general results to several models 
with and without electroweak singlet neutrinos. Concluding comments are given in section 
5. 

2 Theorem on the Number of Physical Phases 

The leptonic mass terms are taken to arise from interactions which are invariant under the 
standard model gauge group Gsm = SU(3) x SU(2) x U(l), via the spontaneous symmetry 
breaking of Gsm- in the standard model and its supersymmetric extensions, the resul- 
tant mass terms appear at the electroweak level via (renormalizable, dimension-4) Yukawa 

4 The NOMAD experiment uses a method proposed in Ref. || while the CHORUS experiment combines 
this with a search for the t tracks in an emulsion. 
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couplings 5 ' 6 . It follows that the mass terms for the charged leptons can be written in terms 
of the Gsm lepton fields as 

3 

- C mass/ = [(L jL hM$ £ kR ] + h.c. (2.1) 
j,k=l 

where j, k = 1,2,3 denote generation indices, and Lj L = (^) L > are the left-handed SU(2) 
doublets of leptons (where £i = e, £ 2 = /i, £ 3 = r). The index a = 2 in (Lj L ) 2 = IjL is the 
SU(2) index. Further, i^R are the right-handed SU(2) singlets, and is the charged lepton 
mass matrix. We use the result here from LEP and SLC that there are three generations of 
standard model fermion with associated light neutrinos []]]. 

In contrast to the charged lepton sector, where one at least knows the relevant fields, in 
the neutral lepton sector, a priori, one does not; in addition to the three known left-handed 
I = 1/2, J 3 = 1/2 Weyl neutrino fields Vji, j = 1 — 3, there could be some number n s of 
electroweak-singlet neutrino fields Xj,R-> j = 1,2, ...n s . The general neutrino mass matrix is 
given by 

1 / M (i) M (D) \ / v c \ 

- Cmass,N = ^ L X° J ^ ( M (0))T M (R) J {^J + h.C. (2.2) 

where v h = (z/ e , z/ M , z/ T )^, \ R — (Xi, ■••■> XnJn! and M^ R ' are 3x3 and n s x n s Majorana 
mass matrices, and is a 3-row by n^-column Dirac mass matrix.^ We denote the 

(3 + n s ) x (3 + ri 8 ) neutrino mass matrix in ( |2.2| ) as M^ N \ where denotes "neutral", 
Q = 0. (In Ref. [0] this was labelled equivalently as M^°\) Thus, 

Mf = Mjp (2.3) 



5 We do not consider models in which fermion masses arise via multifcrmion operators, which arc not 
perturbatively renormalizable. Note that at a nonperturbative level, lattice studies j^] show that a lattice 
theory with a specific multifermion action and no scalar fields may yield the same continuum limit as a 
theory with a Yukawa interaction. 

6 In a number of interesting models, some of these Yukawa couplings are viewed as originating at a higher 
mass scale, such as that of a hypothetical supersymmetric grand unified theory (SUSY GUT) or some other 
(supersymmetric) theory resulting from the E « (a/) -1 / 2 limit of a string theory. It is possible that 
interactions which appear to be Yukawa couplings at a given mass scale, could be effective, in the sense that 
some of the elements of the associated Yukawa matrices could actually arise from higher-dimension operators 
at a yet higher mass scale. The plausibility of such higher-dimension operators can be inferred either from 
the nonrenormalizability of supergravity or as a consequence of the E « (a') -1 / 2 limit of a string theory. 

7 In the context of supersymmetric extensions of the standard model, we assume unbroken R parity so 
that the neutrinos do not mix with the neutralinos (higgsinos and neutral color-singlet gauginos). 
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where j, k = 1,2, 3; 



with j, k = 1, ...n s ; and 



= M%l )k+3 (2.4) 



<> = <ls (2-5) 



where j = 1, 2, 3, = 1, .., n s . In general, all of these matrices are complex. Recall that the 
anticommutativity of fermion fields in the path integral and the property C T = —C (where 
C is the Dirac charge conjugation matrix) imply that 

M (f) = M (f)T ^ j Qr f = l,R (2.6) 

so that 

M (N) = M (N)T (2J) 

Thus, the situation in the leptonic sector is qualitatively more complicated than that in the 
quark sector because of the general presence of three types of neutral fermion bilinears, Dirac, 
left-handed Majorana, and right-handed Majorana, each with its own gauge and rephasing 
properties. The diagonalization of M™' yields the three charged lepton mass eigenstates 
e m , /i m , and r m , while the diagonalization of yields, in general, 3 + n s nondegenerate 

Majorana neutrino mass eigenstates (some could, of course, be degenerate in magnitude, 
leading to Dirac neutrino states). 

A comment is in order concerning possible electroweak-singlet neutrinos. In the super- 
symmetric generalization of the standard model, which is motivated by its ability to main- 
tain the condition v « Mp beyond tree level without fine-tuning, the component fields 
Xj,L wou ld arise as part of Gs^f-singlet chiral superfields (all chiral superfields will be 
written as left-handed). One may recall that, in general, chiral superfields which are singlets 
under the standard model gauge group Gsm can destabilize the hierarchy flQfl . However, in 



commonly used models, the symmetry (such as matter parity) which prevents excessively 
rapid proton decay also excludes the types of terms which would destabilize the hierarchy. 

To count the number of unremovable, and hence physically meaningful, phases in 
and one rephases the lepton fields in ( |2.1|) and (|2.2|) so as to remove all possible phases 
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in these matrices. We gave a general theorem on this counting problem in Ref . J7| . To make 
our discussion reasonably self-contained, we briefly review it here. The rephasings are 



L jL = e~^L' jL (2.8) 
£ jR = & £' jR (2.9) 



for j = 1, 2, 3, and, if there exist any XjR s , & ^ so 



XjR = e^ X) X ' lR (2.10) 

for j — 1, ...n s . In terms of the primed (rephased) lepton fields, the mass matrices have 
elements 

Mf = e ^M £) )M$ (2.11) 
for the charged leptons, and, for the neutrino sector, 

M^' = e^+^M^ (2.12) 



for j,k = 1,2,3; 



for j, k = 1, n s , and 



M^' = e % ^ +p * >M£> (2.13) 

M (DV = e iia j+ ^) M (D) (214) 



for j = 1, 2, 3 and k — 1, n s . 

Thus, if MW has N e nonzero (and, in general, complex) elements, then the Ng equations 
for making these elements real are 

a 3 + = -arg(M^) + (2.15) 

where the set {jk} runs over each of these nonzero elements, and rf^ = or 1.Q 

8 The rjjk term allows for the possibility of making the rephased element real and negative; this will not 
affect the counting of unremovable phases. 
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Similarly, in the neutrino sector, if M^ D \ M^ L \ and M^ R > have, respectively, No, Nl, 
and Nr nonzero (and, in general, complex) elements, then the corresponding equations for 
making these elements real are 

a, + 4* } = -arg(M^) + ffi* (2.16) 

a j + a k = -arg(M^)+ V ^\ (2.17) 

and 

+ 4 X) = -arg(M^) + (2.18) 

where the various ranges of indices are obvious from (2.2). Let us define the (6 + n s )- 
dimensional vector of fermion field phases 

v = ({«,}, {A W },{A (x) }f (2-19) 

where {a,} = {ai,a 2 ,a 3 }, = {P? , , A% and {(3^} = {(3{ x \ 0$}. We also 

define the vector of phases of elements of the various mass matrices^ 

w = ({- axg(iW$) + r^vr}, {{- arg(MW) + £>7r}) T (2.20) 

The dimension of the vector w is equal to the number of rephasing equations N eq = N#+Nn, 
where = + Nr + Np. The indices jk and mn in ( |2.20 ) run over the nonzero elements 



of M™> and the independent nonzero elements of the symmetric matrix M^ N \ We can then 
write (P7HSD-(png|) as 

Tv = w (2.21) 
which defines the matrix T. Since T is an N eq -row by (6 + ra s )-column matrix, clearly 

rank{T) < min{A^ eg , 6 + n s } (2.22) 

Recall that in the quark case M, the analogous inequality was rank(T) < min{N eq , 9}, 
but for realistic models where N eq > 9, one had the stronger inequality that rank(T) < 8 
because one overall rephasing left the Yukawa (or equivalently, mass) matrices invariant. In 



In eq. ( 2.2d|) we change notation slightly from Ref. M, where an overall minus sign was absorbed in T. 



contrast, in the corresponding present case, where N eq > 6 + n s , there is no overall rephasing 
(unless M^ L ) and both vanish) which leaves and M^ N ' invariant, and hence the 

inequality may be realized as an equality, i.e., rank{T) may be equal to 6 + n s . This will be 
illustrated in specific models to be discussed in section 4. 

The first main theorem is then 0: The number of unremovable and hence physically 
meaningful phases N p in and is 

N p = N eq - rank{T) (2.23) 

For the proof, see Ref. J7[. As noted there, this is independent of whether or not the 
matrices are initially (complex) symmetric. There is a one-to-one correspondence between 
each such unremovable phase and the phase of a rephasing-invariant product of elements of 
mass matrices. We discuss these invariants next. 

3 Rephasing Invariants and Theorems on Locations of 
Phases 

A fundamental question concerns which elements of M^' and M^ N ' can be made real by 
rephasings of lepton fields. In Ref. several theorems were given which provide a general 
answer to these questions. The general method is to construct all independent complex 
products of elements of the , / = £, N, having the property that these products are 
invariant under the rephasings (|2.8| )- (|2.10|) . We refer to these as complex invariant products. 
Since in general, by construction, these invariant products are complex, each one implies a 
constraint, which is that the set of elements which comprise it cannot all be made real by the 
rephasings ( ^.8[ )-( PTTD[ ) of the lepton fields. We define an irreducible complex invariant to be 



one which cannot be factorized purely into products of lower-order complex invariants. Since 
reducible complex invariants do not yield any new phase constraints, it suffices to consider 
only irreducible complex invariants, and we shall do so. We define a set of independent 
(irreducible) complex invariants to be a set of complex (irreducible) invariants with the 
property that no invariant in the set is equal to (i) the complex conjugate of another invariant 
in the set or (ii) another element in the set with its indices permuted. This does not imply 
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that the arguments of a set of independent (irreducible) complex invariants are linearly 
independent. We define Ni a to denote the number of linearly independent arguments among 
the independent complex invariants. Further, we define Ni nv to be the total number of 
independent (irreducible) complex invariants in a given model. It is useful to define Ni nVi 2 n 
to denote the number of independent complex invariants of order 2n. 

In Ref. |7J we discussed the general construction of complex rephasing invariants, enu- 
merated the various types of these invariants at fourth and sixth order, and gave explicit 
expressions for the fourth order and a number of sixth order complex invariants. Here we 
shall give a general set of explicit expressions for complex invariants up to sixth order, inclu- 
sive. We recall our classification of rephasing invariants into two general types: (1) P-type, 
which involve elements of mass matrices in one charge sector (charged leptons or neutrinos), 
and (2) Q-type, which involve elements of mass matrices in both of these charge sectors. We 
recall that the two quartic P-type invariants are 

pif) _ pif) _ M U) M if) M iJ) *M iJ) * (1 1) 

' k n k :J .h,y: ~ r jlki,j 2 h 2 ~ IV1 jik 1 IV1 j 2 k 2 IV1 j2k 1 IV1 j 1 k 2 i 6 - 1 ) 

for / = £,N. Note the symmetries 

pif) _ pif) (a 2) 

pif) = pif)* to o\ 

j\k\hk 2 jik2,j2ki \°-°) 

(whence also PjHijafa = Pjlkthki) ^ or / = ^ ^ anc ^ ^ or / = ^ ne additional symmetry 

P {N) • • = P (iV) ■ • (3.4) 

0132,3334, 3231,3433 V. u *^7 

In a convenient shorthand notation, the complex invariants in ( |3.1| ) for f = £,N may be de- 
noted as £££*£* and NNN*N*, respectively. The P 4 complex invariants of the form NNN*N* 
may be further classified into six subtypes according to which submatrices in they in- 

volve, M^ L \ M( r \ and/or M^: 

NNN*N* : LLL*L* LDL*D* 
RRR*R* RDR*D* 

DDD*D* LRD*D* (3.5) 



(It is easily seen that the rephasing invariants of the subtype LRL*R* are real.) Explicitly, 
the first three complex invariants are given by (|3.1|) for / = L,R,D, respectively, and the 
other three by the products 

Tl(LD) _ M (L) M (D) M (L)* M (D)* fi \ 
n jij2,j 3 ki - 1V1 jij2 1VI j3ki m jij3 1V1 j2k! V 6 - ) 



agSw. = <>2<»2' (3-7) 

and 

q(LRDD) _ M (L) M (R) ,,(£>)* M (D)* ,„ ^ 

It is straightforward to express these in terms of the elements of P( N \ using (|2.3|) - ([2.5|) : 

p (L) p {N) ,„ q) 

3132,3334. ± 3132,3334. V°- u I 



Pk^Mki — k 2 +3, k 3 +3 k 4 +3 (3.10) 



p{D) _ p {N) , , 

r 3iki,hk 2 — j 2 k 2 +3 \3- Ll ) 

n (LD) = p (N) ,„ m 
3l32,33ki 32 31, 33 kx+3 ) 

W (RD) _ p (N) /o I o\ 

^kik 2 ,j!k 3 ~ ^fei+3 k 2 +3, ji k s +3 {O.LO) 

a (LRDD) _ p(AT) /o 

lL h32,kik 2 — r j 2 h, ki+3 k 2 +3 K^-^J 

At quartic order, we found two Q-type complex invariants, namely 

m _ M (D) M {i) M (D)* M {i)* f o 1 r\ 

^jiki^rri! — IVI jik 1 IVI j 2 m 1 IVI j 2 k 1 IV1 jimi I^.IDJ 

and 

3132,33^1 J1J2 33™l JU3 J2ini ^ u u ^ 
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In the short notation, these are denoted D£D*£* and L£L*£*. Note that 



tt( l /) _ f)(^'J > 

3132,33ml ^hhJami 



(LJ) 



(3.17) 



At sixth order there are first the complex invariants each of which involves only one 
charge sector, ////*/*/* for f — £,N. The explicit expressions for these are 



P 



M U l M U) h M if l M U) h * M U) * M U l* 



(3.18) 



for f = £,N. For any /, the invariants in (3.18) satisfy 



P U) 
jifti Jatea >J3&3 



p U) 



p U) 

33^3 Jlfcl J2fc 2 



and 



p (f) = p(f)* 

jiki jate ijafa ^3^2 J2fci ,ii ^3 



(3.19) 



(3.20) 



For the cases where is symmetric, namely Af^ = Af'^' or, for submatrices, M^ L \ M^ R \ 
one has the additional symmetry 



P (f) P (f) 

hh ,3304, ,3536 32 jl ,3635 J4J3 



(3.21) 



The NNNN*N*N* invariants may be divided into 11 subtypes according to which subma- 
trices of M^ N > they involve: 



NNNN* N*N* 



LLLL*L*L* 

D D D D* D* P* 
tttttttt it it 

DDDD*D*D* 
DDLD*D*L* 
DDRD*D*R* 



LLDL* L* D* 

RRDR*R*D* 

LRDL*R*D* 

LLRL* D* D* 

LRRD*D*R* 

LRDD*D*D* 



(3.22) 



Complex conjugates of these invariants yield the same phase constraints and hence are not 
listed. Further, in the shorthand notation, the ordering of the factors is not important; 
e.g. DDLD*D*L*, DLDD*L*D*, and LDDL*D*D* represent the same subtype, and so 
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forth for the others. Note that the 6'th order rephasing invariants £RR£*R*R*, ££R£*£*R*, 
LRRL*R*R*, LLRL* L* R* , and £LR£*L*R* all reduce to real factors times quartic complex 
invariants and hence are not irreducible 6'th order complex invariants. Given an explicit 
formula for a particular type of NNNN*N*N* complex invariant in terms of submatrices of 
, it is straightforward to use eqs. (|2.3| )- (|2.5| ) to reexpress it in terms of M^ N \ Thus, for 



LLLL* L* L* , we have the relation P^ hM = P^Lu,^ DDDD*D*D* the relation 

P hkihk 2 , j3 k 3 = ^2+3,^+3,^3+3' alld SO forth fOT the ° therS - 

Explicit expressions for the other subtypes of NNNN*N*N* complex invariants in ( |3.22| ) 
are given below (in some cases there is more than one kind of structure for a given subtype): 

DDLD*D*L* : M^M^M^M^* 'M%£ 'M%£ (3.23) 



M {D) M (D) M (L) M (D)* w(fl> M (L)* /„ 



DDRD*D*R* : M^2<IM2 ^ M Sl ^ 

^M^Mj^Mj^M^M^ (3.26) 

LLDL* L* D* : M^M^M^M^M^M^ (3.27) 

RRDR*R*D* : ^CCCXtC ( 3 ' 28 ) 

LRDL*R*D* : M>SX£C (3.29) 

LLRL* D* D* : M^M^M^M^ (3.30) 

LRRD* D* R* : M^MfflXlK£ Kk\* ( 3 - 31 ) 
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and 



LRDD* D* D* : M^M^M^Mj^ M%> M%> (3.32) 

M (L) M (R) M (D) M (D)* M (D)* M (D)* /„ oo\ 

Note that the invariants in (|3.23|) and ([3.27 ) are of the form 



nVff) _ M U) M (f) M (f) M (f) *M (f) *M (f,) * n 34) 

for (///0 = (DDL) and (LLD) , respectively. 

We also find the following seven independent types of 6'th order complex invariants 
linking the charge Q = — 1 and Q = sectors: 

ULttL* IDL£*D*L* 

LL£L*L*£* £DR£*D*R* 

UDttD* £LR£*D*D* 
DD£D*D*£* (3.35) 



The explicit forms for these are 



M W M {!) M (L) M (i)* M (i)* M (L)* ( ~ o fi \ 

Jimi J2»2 J3J4 Jim2 J3"ll J234 yu.vvj 



M (t) M (n M (L) M (i)* M (i)* M (L)* ( ~ o 7 \ 

1VI j 1 m 1 1V1 j 2 m 2 m hh m hmi 1V1 j4m 2 IVI jij2 l " ' ) 



LL£L*L*t : ^S^CCCC ( 3 - 38 ) 



«I>rrD- : Mj^J^^^Mj^M^ (3.39) 



£DR£*D*R* : ^Mg^^^M^* ( 3 - 40 ) 
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£LR£*D*D* 



(3.41) 



£DL£*D*L* : MJKi^M^CC ( 3 - 42 ) 



M (£) M (D) M (L) M {1)* M (D)* M (L)* , 
M jim 1 1VI j2k 1 1VI j3H IV1 j4m 1 1V1 j 1 k 1 M hh \6A6) 



Mf mx Mfl MWM^ M^i (3 .44) 



M (i) M (D) M (L) M (£)* M (D)* M (L)* (14.^) 



M (() M (D) M (L) M {1)* M (D)* M (L)* ,„ 4fi \ 



M (i) M (D) M (L) M (i)* M (D)* M (L)* 

IVI jxm 1 IVI j 2 k l IVI j3M IVI j4mi IVI j 3 ki m jij 2 W-^' ) 

The methods for determining these are similar to those discussed in detail for the quark case 
in Ref. 0. For example, the general explicit form of the £DL£*D*L* invariants is 

where a G S4, the group of permutations of 4 indices. We then determine the full subset of 
these 4! products which yields independent irreducible sixth order complex invariants, and 
obtain the results listed in (p^ -([347D. Note that (|3~36|) and (|3~38| ) are again of the form 
( Egg ) for (///') = (££L) and (LL£), respectively. 



In the case of quark mass matrices (with Nq = 3 generations of quarks), we showed || 
that the fourth and sixth order complex invariants sufficed to describe all of the phase 
constraints on the mass matrices. An analogous result does not hold in the leptonic sector, 
however; in general, it may be necessary to include complex invariants of higher order than 
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sixth to obtain all of the phase constraints. The reason for this is just that the neutral part 
of the lepton sector is, in general, more complicated than either the up or down quark sector. 
Rather than enumerate in general all of the possible eighth-order complex invariants, which 
would be quite involved, we will make clear how one uses them in practical calculations via 
two specific models in section 4. 

We next review our theorems on invariants, as applied to lepton mass matrices. For 
a given model, one constructs the maximal set of N ia independent complex invariants of 
lowest order(s), whose arguments (phases) are linearly independent. Then (a) each of these 
invariants implies a constraint that the elements contained within it cannot, in general, all 
be made simultaneously real; (b) this constitutes the complete set of constraints on which 
elements of and can be made simultaneously real; and hence (c) N p = N ia . 

An immediate corollary is: If in a given model there are as at least as many independent 
quartic complex invariants with independent arguments as there are unremovable phases, 
then the arguments of all higher order complex invariants are expressible in terms of those 
of the quartic invariants and hence yield no new phase constraints. 

In the quark case, we presented a useful graphical representation for the complex invari- 
ants in Ref. ||- This again works directly for invariants involving only the charged lepton 
sector. For invariants involving the neutral lepton sector, there is no 1 — 1 correspondence 
between a matrix element M-^ occuring in a given complex invariant and a point on an 
graphical array representing the matrix M^; instead, as a result of the symmetry (|2.7|) , 
M(n)t _ m(N) ^ there is a obvious 2- fold homomorphism according to which the element 
= Mf,j corresponds to the two points (j'th row, /c'th column) and (A;'th row, j'th 
column) in the array representing M^ N >. For an invariant involving 2n elements of M^ N ' of 
which a (possibly null) subset of r are diagonal elements, there are 2 2n_fc -ways of representing 
it graphically. Among these ways, however, one can always find a graphical representation 
which is analogous to those for the quark sector complex invariants discussed in Ref. [Q. 
Moreover, as a tool for quickly determining the complex invariants, given the forms of the 
mass matrices, the graphical method is equally useful in the leptonic case as in the quark 
case. We shall illustrate it in this context below. 
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Since the full set of N{ nv independent complex invariants will have arguments which are 
not, in general, linearly independent, it follows that 

N inv > N ia (3.49) 

It may also happen that, e.g. for order 2n = 4, the number of independent quartic complex 
invariants, N inVti is greater than N ia . Using the methods of Ref. (see also Ref. 0), one 
can determine how many of the complex invariants of all orders, or of a particular order, 
have linearly independent arguments, and can thus construct sets of such invariants with 
this property. 

4 Applications to Specific Models 

Although our theorems are quite general, it is useful to see how they apply to various 
specific models. In Ref. we gave a brief discussion of applications. Here we will consider 
several more models. In contrast to the quark sector, where one has reasonably accurate 
measurements of most Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix elements, 
in the leptonic sector, there is not even any definite direct evidence of leptonic mixing, let 
alone accurate measurements of various mixing matrix elements. As noted above, probably 
the strongest indirect evidence comes from the apparent deficiency of the solar neutrino flux. 
The situation concerning possible evidence for atmospheric neutrino oscillations is unclear 
at present [p], ||. We shall tentatively assume that the solar neutrino deficit does indicate 
neutrino oscillations.^ Moreover, since our purpose is to illustrate the application of our 
general results on phases, we will not discuss the details of the phenomenology of the models. 

In general, models for fermion masses and mixing can be classified according to whether 
they assume a theoretical framework of perturbative or nonperturbative electroweak sym- 
metry breaking. We shall concentrate here on the class of models in which the observed 
electroweak symmetry breaking is perturbative. In this class, an appealing framework is pro- 
vided by supersymmetric extensions of the standard model, which stabilize the Higgs sector 

10 For an recent phenomenological study which also includes fits to atmospheric neutrino oscillations, see, 
e.g., Ref. 0. 
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and hence stabilize the hierarchy between the electroweak scale and the Planck scale. ||12 
In these models, one typically hypothesizes some simple forms at a high mass scale (e.g., a 
mass scale characterizing a possible grand unified theory (GUT) or a mass scale near to the 
string scale, Mp) for the Yukawa and higher- dimensional operators which are responsible for 
fermion mass generation. As part of this, one assumes some symmetries to prevent various 
Yukawa couplings and higher-dimensional operators, and thereby render various entries in 
the (effective) Yukawa matrices zero. The purpose of this is, of course, to minimize the num- 
ber of parameters and hence increase the apparent predictiveness. One then evolves these 
forms for the Yukawa matrices down to the electroweak mass scale using the appropriate 
renormalization group equations. For this evolution, we use the renormalization group equa- 
tions of the minimal supersymmetric standard model (MSSM). Parenthetically, we recall 
that when one evolves these forms for the Yukawa matrices down to the electroweak level, 
the zero entries do not, in general, remain zero. Clearly, the entries which are modelled as 
being exactly zero might well be nonzero, but small quantities. For example, one possibility 
is that they are suppressed by certain powers of (M r /Mp), where M r is the reference scale 
at which one analyzes the effective Yukawa interaction. Since our results do not depend on 
whether or not the mass matrices are symmetric at some mass scale, we shall consider the 
general case of non-symmetric mass matrices. 

4.1 Model Without Electroweak-Singlet Neutrinos 

We first consider a model with no electroweak-singlet neutrinos, i.e., with n s = 0. This is 
defined by [§. 
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E23 







E32 


E33 
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L>12 


L22 


L23 







L23 


L33 



(4.1.1) 



= L 12 L22 L23 (4.1.2) 
\ L23 L33 J 

(where each element is, in general, complex, since no symmetry forces it to be real). We 
have verified that, for appropriate choices of the parameters, this model is able to fit current 
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solar neutrino data and other established limits on neutrino masses and mixing. For this 
model, there are N eq = 10 complex elements in and and corresponding rephasing 

equations. We calculate the matrix T and find rank(T) = 6. Our theorem ( |2.23j ) then 
implies that there are N p = N eq — rank(T) = 4 unremovable phases in M^' and M^ L \ 
(This model provides an illustration of how rank(T) may realize the inequality (|2.22|) as an 
equality.) We find N inv ^ = 8 independent complex (quartic) invariants: 

^2,33 = -^22-^33-^32-^23 (4.1.3) 

-^22,33 = -^22-^33-^23 (4-1-4) 
1412,22 = -^12-^22-^22-^12 (4.1.5) 
1^12^32 = -^12-^32-^23-^12 (4.1.6) 



n22^32 — -^22-^32-^23-^22 (4-1-7) 



1^22^33 = -^22-^33-^23-^23 (4.1.8) 
= -^23-^32-^33-^22 (4.1.9) 
^23% = -^23 £33 £33-^23 (4.1.10) 

From these we calculate the corresponding Z matrix and find that it has rank 4, so that of 
the eight arguments of the complex invariants, there are A^ Q = 4 linearly independent ones, 
in accord with the equality iVj a = N p and our result that N p = 4. There are correspondingly 
four independent phase constraints. In particular, for this model it is not, in general, possible 
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12,22) ii 22,32J 



by any rephasings to make either the charged lepton or neutrino mass matrices real. A set 
of complex invariants with linearly independent arguments is given by 

(4.1.11) 

We denote the arguments of the invariants in eq. ( |4.1.11| as 8j, j = 1,...,4; i.e., 9\ = 
arg(P22ss), etc. The only element which does not appear in any complex invariant and 
hence may be rephased freely is E^i- 

A necessary task which a model builder must carry out is to determine which of the 
elements of the mass matrices can be rendered real by the rephasing of the lepton fields. We 
display below an example of an allowed form for the mass matrices after such rephasing is 





I -^21 1 
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1-^12 I 


L22 
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This model has iV, 



invfi 



4 sixth order complex invariants: 

Q22,3342 — -^22-^33-^12-^23-^12-^32 



(4.1.12) 



(4.1.13) 



(4.1.14) 



n {UL) 
V22,33,12 



n (LL£) 
V33, 12,22 



MIL) 



771 771 r 771* 771* T* 
-£^22-^33-^12-^23-^12-^23 



T T Z? T *2 TP* 
-^33-^12-^22-^23-^12 



771 771 r 77* 77* T* 

-^33-^12-^22-^32-^23-^12 



(4.1.15) 
(4.1.16) 
(4.1.17) 



^33,12,22 

However, as a result of our corollary 1 above, since we have already exhibited a set of quartic 
complex invariants with independent arguments and comprised of N ia members equal to 
the number N p of unremovable phases in the mass matrices, it follows that all higher order 
invariants, and, in particular, the sixth order ones listed above, have arguments which can 
be expressed as linear combinations of those of the set of N ia fourth order invariants and 
hence do not yield any new phase constraints. 



4.2 Models with Electroweak- Singlet Neutrinos 
4.2.1 n s = 3 Model 1 

We will concentrate on models with n s = N G = 3. Our methods are easily applied to models 
with other values of n s . We begin with an idealized model in which left-handed Majorana 
mass terms are assumed to be negligible, i.e., = 0, and the other mass matrices are 

given by 
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^33 , 





We have again checked that with appropriate choices of the parameters, this model can 
fit the experimental constraints discussed at the beginning of the paper. The model has 
N eq = 15, and we compute rank(T) = 9, so that there are N p — 6 unremovable phases in 
the mass matrices. The model has N inv ^ = 6 quartic complex invariants: 

^22% — -^22 £33 E23 El 2 (4.2.4) 



Qt3,% — -^23-^33-033-^23 



Q?3% — -D23-E32-D33-E22 



Qv2% — -D12-E32-D32-E12 



p n z?* n* 

"23,33 _ -^23 -^33-^33 -^32 
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(4.2.5) 



(4.2.6) 



(4.2.7) 



(4.2.8) 



and 

^22^33 = -^22-^33-^23 (4.2.9) 

Using our Z matrix method, we find that among these six invariants there are five linearly 
independent arguments. A set of five quartic complex invariants with linearly independent 
arguments is provided by two among the subset {P 22j33 , $23,33) Q^%} together with the 
three invariants Qi2%, £23,33, anc ^ ^22^33 ■ To obtain an invariant encoding the sixth unre- 
movable phase, we consider 6'th order complex products. We find N inVi6 = 7 of these: 

I 6a = E 12 E 23 D 32 E* 2 E* 33 D* 12 (4.2.10) 



I 6b = D 23 D 32 E 12 D* 12 D* 33 E* 22 (4.2.11) 

he = E 12 D 23 R 22 E 22 Dl 2 Rl 3 (4.2.12) 

hd = E 12 D 33 R 22 E 32 Dl 2 R 23 (4.2.13) 

he = E 12 D 23 R 23 E* 22 D* 12 R* 33 (4.2.14) 



hf — E 12 D 33 R 23 El 2 D* l2 R* 33 (4.2.15) 

and 

I Qg = D 2 23 R n D*lR* 33 (4.2.16) 

Among these seven 6'th order invariants, the arguments of the first six are linearly dependent 
upon the five arguments from the quartic invariants, but the argument of h g is linearly 
independent of these and hence accounts for the sixth unremovable phase in the model. 
Therefore a complete set of complex invariants is given by 

snW n( D V ~ r 1 (a 9 1 7^ 

lV23,33) V23,32> V 12,32) "23,33 5 - r 22,33) 1 &gi yt.A.i-i) 
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We denote the arguments of these complex invariants as j = 1, 6, i.e., X = arg^^-^), 
etc. An example of an allowed form for the mass matrices in this model after rephasing of 
the lepton fields is 



# = 



M (RY 




(4.2.18) 



(4.2.19) 



and 
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I d 2 i 1 




\Dn\ 


\D 32 \ 



(4.2.20) 



\ 

\D 23 \ 
I-D33I / 

By considering a special case of this model, we can illustrate our statement in section 
3 that for leptonic mass matrices, in contrast to quark mass matrices, the quartic and 
sixth order complex invariants do not necessarily suffice to encode all of the unremovable 



phases. For this purpose, let us, for example, set -R33 = 0. In this case, N t 



eq 



14 and 



the rank of the resultant T is again 9, so that N p = 5. Among the quartic invariants, the 
first four remain, but the last two, 2^33 and P22I3, evidently vanish. From our previous 
analysis it follows that among the remaining four quartic invariants, there are three linearly 
independent arguments. Of the original seven sixth-order complex invariants, the last three 
evidently vanish. Since the arguments of the first four sixth order invariants are all linearly 
dependent upon those of the quartic invariants, they does not account for either of the two 
remaining unremovable phases. Consequently, one must proceed to construct the complex 
8'th order complex invariants for this special case. We find the following two: 



ha — E22D23D 32 RiiE^ 2 D2lRl 3 



(4.2.21) 



^86 — -^23-023-032-^11 -^33-021-^23 (4.2.22) 

The arguments of these two invariants are linearly independent of the three arguments from 
the quartic invariants, yielding the full set of N ia = N p = 5 independent arguments for 
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this special case of the model where R33 = 0. A set of complex invariants with linearly 



independent phases is provided by (i) three out of eqs. fl4.2.4D - (|4.2.9T) together with (ii) 
I$ a and Isb- (The general model with R 33 7^ has several additional 8'th order complex 
invariants, but they do not yield any new phase constraints.) 

4.2.2 n s = 3 Model 2 

A second model with n s = 3 electroweak-singlet neutrinos is defined by the mass matrices 

(4.2.23) 



(4.2.24) 
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(4.2.25) 



(4.2.26) 



This model has N eq = 15. We calculate the 15 x 9 matrix T and find rank(T) = 9, so 
that, according to our theorem, there are N p — 6 unremovable phases in the leptonic mass 
matrices. To determine which elements can be made real, we construct a set of complex 
invariants which yields the full set of phase constraints. We find that there are Ni nVt 4 = 8 
quartic complex invariants: 

(4.2.27) 



Vl2,32 



-D12-E32-D32-E12 



^23,33 — -^23 £-33-033-^23 



^(RD) 75 j}Sf p,* 

"21,23 — -"-12 -^23-^23 JJ 21 



1 33^32 — -R33-D32-R23-O33 



(4.2.28) 
(4.2.29) 
(4.2.30) 
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^{RD) d n r?* n* 

"21,32 — -"42 -^32-^23 -^21 



(4.2.31) 



(LRDD) 
12,12 



Li2Ri2D 2 iD 12 



(4.2.32) 



(LRDD) 
12,23 



-^12-^23-023-^12 



(4.2.33) 



n (LRDD) 



6 33,23 



l 33 r 23 d; 2 d; 3 (4.2.34) 

By our usual Z matrix method, we determine that among these eight invariants there are 
six linearly independent arguments, so that the quartic complex invariants suffice to yield all 
phase constraints in this model. We find that a subset of the eight quartic invariants which 
have independent arguments is given, e.g., by 



r n {Dl) n ipi) ^(RD) n (LRDD) n (LRDD) n (LRDD) x 
1V12,32? V23,33? "33,32; "12,12 j "12,23 j "33,23 J 



(4.2.35) 



We denote the arguments of these complex invariants as uuj, j = 1,..,6, respectively. An 
example of an allowed set of mass matrices after rephasing is 
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(4.2.36) 



(4.2.37) 



(4.2.38) 



\D l2 \ \ 
) 2 i|e"^ 5 \D 23 \e- iuj( ' (4.2.39) 
\D 32 \ \D 33 \ J 

By our corollary 1, since we have constructed a set (|4.2.35|) of quartic invariants with the full 

Nia = N p set of linearly independent arguments, all higher order complex invariants have 
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arguments which can be expressed in terms of those of this set, and hence they do not yield 
any new phase constraints. For completeness, we exhibit the 6'th order complex invariants 
in this model. Two of these are of LRDD*D*D* type, 

^12^6,35 — Li 2 R 33 D 32 Dl 2 Dl 3 D 33 (4.2.40) 



^26^33,54 — -^33-^12-023-021-^32-^33 



(4.2.41) 



while the third is of DDLD*D*L* type: 



^12*63,35 — D 32 D 33 Li 2 Dl 2 D2 3 Ll 3 



(4.2.42) 



(Several choices of indices on the Pj^j^jz^ yield the same invariants; we have listed only 
one for each case.) 



4.2.3 n s = 3 Model 3 



We next give two toy models to illustrate some theoretical points. The first is an example 
of a model with no complex quartic invariants, where the first complex invariant occurs at 
sixth order. It is defined by 



(4.2.43) 
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(4.2.44) 



(4.2.45) 



(4.2.46) 



24 



We calculate that rank(T) = 9 for this model, so that there is N p = 1 unremovable phase. 
The single independent complex sixth order invariant is of DDRD*D*R* type: 

D 2 32 R 33 D* 3 lR* 22 (4.2.47) 

The corresponding phase may be placed in any of the mass matrix elements involved in 



4.2.4 n s = 3 Model 4 



This model illustrates that, in contrast to the case in the quark sector, the quartic and 
sixth order complex invariants do not necessarily suffice to yield all phase constraints in the 
leptonic sector. The model is defined by and as in eqs. ([4.2.43 ) and (|4.2.46|) , 

respectively, together with 



M (L) =| L 22 | (4.2 

and 
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R22 I (4.2.49) 




V 

The T matrix for this model has rank 9, so that N p — 1. There are no (nonzero) complex 
quartic or sixth order invariants. The model has the single independent 8'th order invariant 

L 22 R 22 Dl 3 D* 2 \Df 2 (4.2.50) 

As before, the corresponding phase may be placed in any of the mass matrix elements 
involved in (|4.2.47|) . 



4.3 Unified Models of Quark and Lepton Mass Matrices 

One of the appeals of grand unified theories has always been the fact that they allow one 
naturally to relate quark and lepton Yukawa, and hence mass, matrices.^ In this context, we 

11 This appeal is offset somewhat by the fact that even when one considers supersymmetric grand unified 
theories to stablize the overall gauge hierarchy, there is a still a serious problem associated with the necessity 
of splitting the electroweak doublet Higgs fields from color triplet Higgs fields which occur in the same 
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note that the methods for the counting and allowed placements of the unremovable phases in 
the lepton sector discussed here and in Ref. [0 can be combined with our analogous methods 
for the quark sector discussed in Ref. U in the context of these unified models of lepton and 
quark masses. If various elements of the quark and lepton mass matrices are equal, or differ 
only by some real factor (because of the types of Yukawa couplings allowed by the grand 
unified gauge invariance and by various other symmetries which are imposed to restrict these 
couplings), some of the invariants in the quark and lepton sectors may coincide, leading to 
a reduction in the total number of phases. We gave an example of this in Ref. [|7| . 

5 Conclusions 

The goal of understanding fermion masses and mixing remains one of the most important 
outstanding problems in particle physics. In this paper we have given a detailed discussion of 
our methods for determining the number of unremovable phases in a given model of leptonic 
mass matrices, including new results on higher order invariants and illustrative applications 
to specific models. As ongoing and future experiments searching for evidence of neutrino 
masses and lepton mixing yield new information, the construction of more tightly constrained 
predictive models of fermion masses will continue. We believe that the general tools discussed 
here will be of use in these studies. 

This research was supported in part by NSF grant PHY-93-09888. 



representations of the GUT gauge group. Furthermore, the most likely source for the additional symmetries 
which model-builders impose to get zeroes in the (effective) Yukawa matrices is an underlying string theory. 
But, at least at the limited level that the phenomenology of string theories is understood, they do not 
necessarily yield simple gauge groups (the essence of GUT's) or the requisite adjoint representations of the 
Higgs chiral superfields which are usually necessary for the breaking of the GUT. A typical direct product 
gauge group resulting as the low-energy limit of a string theory can avoid the doublet-triplet splitting problem 
in GUT's. Although one can get GUT's and adjoint (and other higher-dimensional) representations from 
strings by using higher-level Kac-Moody algebras on the string worldsheet, this reintroduces the doublet- 
triplet splitting problem. Of course, one still believes that at a deep level, quark and lepton mass matrices 
should be related in a simple way. 
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